Frequency metrology of helium around 1083 nm and determination of the nuclear 

charge radius 
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We measure the absolute frequency of seven out of the nine allowed transitions between the 2 3 S and 
2 3 P hyperfine manifolds in a metastable 3 He beam by using an optical frequency comb synthesizer- 
assisted spectrometer. The relative uncertainty of our measurements ranges from 1 x 1CP 11 to 
5 x 10 -12 , which is, to our knowledge, the most precise result for any optical 3 He transition to date. 
The resulting 2 3 P-2 3 S' centroid frequency is 276 702 827 204.8 (2.4) kHz. Comparing this value with 
the known result for the 4 He centroid and performing ab initio QED calculations of the 4 He- 3 He 
isotope shift, we extract the difference of the squared nuclear charge radii Sr 2 of 3 He and 4 He. Our 
result for 5r 2 — 1.074(3) fm 2 disagrees by about 4<r with the recent determination [R. van Rooij et 
al, Science 333, 196 (2011)]. 



Spectacular progress of experimental techniques, 
achieved in the last decades, has made precision spec- 
troscopy of light atoms a unique tool for the determina- 
tion of fundamental constants and properties of atomic 
nuclei. The underlying theory, quantum electrodynam- 
ics (QED), allows one to calculate atomic properties ab 
initio and keep control of the magnitude of uncalculated 
effects. Possible discrepancies between theory and ex- 
periment may signal a lack of our knowledge of details 
of the interactions between electrons, nuclei, and other 
particles. An important recent example is the discrep- 
ancy of the proton charge radius derived from the spec- 
troscopy of the electronic and muonic hydrogen [I] . This 
discrepancy is still unresolved and might lead to impor- 
tant consequences, such as a change of the accepted value 
for the Rydberg constant (which was, up to now, consid- 
ered to be one of the best known fundamental constants) 
or discovery of unknown effects in the electromagnetic 
lepton-nucleus interaction. 

Another important disagreement reported recently 
concerns the charge radii of helium isotopes. Specifically, 
the difference of the squares of the charge radii of the 3 He 
and 4 He nuclei determined from the 2 1 £-2 3 5' transition [2} 
was reported to differ by about 4 standard deviations (a) 
from that derived using the 2 3 P -2 3 S transition [3j. In 
this Letter, we aim to resolve this discrepancy by per- 
forming an independent, high-precision measurement of 
the 4 He- 3 He isotope shift, on the one hand, and by ad- 
vancing the theory of the helium isotope shift, on the 
other hand. 

In the experimental part, we measure the absolute fre- 
quency of seven out of the nine allowed transitions be- 
tween the 2 3 S and 2 3 P hyperfine manifolds of 3 He with 
a precision ranging from 1 x 1CP 11 to 5 x 10~ 12 . To the 
best of our knowledge, these are currently the most ac- 
curate measurements for any optical 3 He transition. In 



the theoretical part, we perform a rigorous QED calcu- 
lation of the isotope shift of the centroid of the 2 3 P-2 3 S 
transition, identifying several corrections omitted in the 
previous studies and carefully examining the uncertainty 
due to uncalculated effects. Combining the experimental 
and theoretical results, we obtain the difference of the 
squares of the charge radii of the 3 He and 4 He nuclei. 
The improved theory is also applied for reexamination of 
the previous experimental result [2j. 

The structure of the 2 3 S and 2 3 P levels in 3 He and 4 He 
is shown in Fig.[TJ In the following, we will use the short- 
hand notation P^ s p Fp to denote the 2 3 Si t F s — ¥2 3 Pj Pi f p 
hyperfine transition in 3 He and Pj p to denote the 
2 3 Si — >2 3 Pj p transition in 4 He. We now address the 
measurement procedure applied to the seven transitions 
between the 2 3 S and 2 3 P hyperfine manifolds of 3 He. 

1/2 T-»3 /2 

The remaining two allowed transitions, P 1 3 / 2 anc ^ ^2 3/2' 
have a very weak intensity due to the hyperfine suppres- 
sion 0, [B[ , which prohibits their measurement with our 
spectroscopic setup. 

Multircsonant precision spectroscopy was performed 
by using the optical frequency comb synthesizer (OFCS)- 
assisted laser system at 1083 nm described in Ref. [6j. 
In each run, quasisimultaneous saturation spectra of two 
out of these seven hyperfine transitions were recorded 
in an absolute frequency scale, and the frequency center 
was measured by a fitting procedure [f|. In this way, 
absolute frequencies of each transition and their differ- 
ences are measured minimizing possible time-dependent 
systematic errors. This procedure is repeated for about 
200 runs for each transition. During these runs, differ- 
ent transitions were coupled to each other, and for each 
couple the multircsonant probe laser system was inter- 
changed, in order to randomize, as much as possible, the 
measurement procedure. The final results, together with 
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FIG. 1. (color online) Level scheme of the 2 3 S and 2 3 P man- 
ifolds of 3 He and 4 He isotopes. 3 He HFS, 4 He fine structure 
(FS) and 4 He - 3 He IS splittings are shown. 



TABLE I. Absolute frequency measurements of Pj^ F 3 He 
hyperfine transitions: statistical results and systematic error 
budget, in kHz. Uncertainties are given in parentheses. 



Transition 



Statistical valuc a 



Zccman 
shift " 



FinaP 
result 



1,1/2 
D 3/2 

2,5/2 
D l/2 

2,3/2 
□3/2 

1,3/2 
□3/2 

1,1/2 
□ 1/2 

0,1/2 
□3/2 



276698 164610.4(1.8) (0.233) 276 698164 610.4(2.0) 

276698611209.1(3.1) (0.148) 276 698 611209.1(3.2) 

276 698 832 617.9(2.5) (0.303) 276 698 832 617.9(2.5) 

276 700 392 099.8(0.9) (0.369) 27 6 700 392 099.8(1.3) 

276 704 904 311.7(2.3) (0.233) 276 704 904 311.7(2.4) 

276 726 257468.9(1.1) (0.932) 276 726 257468.9(1.7) 

276 732 997170.4(2.3) (0.466) 276 732 997170.4(2.5) 



a Each measurement was corrected by the day-by-day 

RS+LS+2DS shifts (see text for details). 
b Because of residual magnetic fields (< 0.03 £tT). 
c The OFCS error of 10~ 12 xstatistical value (kHz) and the IDS 

error of 0.8 kHz were added in quadrature in the final 

uncertainty. 



the error budget, are summarized in Table [IJ 

Particular attention was paid to single out and quan- 
tify possible systematic errors. As in our previous 4 He 
measurements 0, Q , the main systematic correction was 
due to recoil- induced mechanical shift (RS) Q. As in 
that case, we calculate its contribution for each transi- 
tion by solving the atomic Bloch equations, including ac 
Stark shift [light shift LS], and taking into account the 
atom dynamics during interaction with the laser, in the 
present experimental conditions. Second-order Doppler 
(2DS) correction due to the longitudinal velocity distri- 
bution of the atoms in the 3 He metastable beam was 
also included in this calculation. Since RS is like an ac- 
cumulated recoil during laser-atom interaction time, it is 
strongly dependent on the metastable atomic flux in our 



experimental setup, and hence on the dc discharge con- 
ditions used to metastabilize the 3 He atomic beam. We 
have noticed that such conditions changed during a day 
of measurements, due to progressive saturation of the fil- 
tering system for contaminant gases, used in the 3 He gas 
recycling line inserted in the atomic beam. In fact, we 
monitored this change by measuring the atomic longitu- 
dinal velocity distribution behavior during a day. From 
this data, an averaged longitudinal velocity distribution 
for each day is determined, which enters as a parame- 
ter in the RS+LS+2DS calculation. All measurements 
for each transition are corrected by the corresponding 
day shift. The final frequency is calculated as a statisti- 
cal average of all corrected measurements performed for 

each transition. Such a procedure is shown in Fig. [2] 
1 li 

for the -P 2 3/2 transition, where about 180 measurements 
without (squares) and with (circles) RS+LS+2DS day- 
shift corrections are reported. As a result, a Gaussian 
distribution of the corrected measurements is shown in 
the fill-bar graph of Fig|3J witnessing that our data have 
"white" statistical fluctuations. 

Systematics uncertainties due to first-order Doppler 
shift (IDS), OFCS accuracy, and Zeeman shift, have been 
added in quadrature to the statistical one, as shown in 
Table HI IDS was avoided due to the saturation spec- 
troscopy configuration but with an error of 0.8 kHz, 
due to the achievable angular accuracy between the for- 
ward and the backward interacting laser beams. The 
Rb-GPS disciplined quartz oscillator, used in our OFCS, 
guarantees a relative accuracy of 10 -12 , considered in the 
error budget. Finally, a residual magnetic field in the 
atom-laser interaction region lower than 0.03 /iT, gives 
a Zeeman shift uncertainty for each transition (Table HI 
third column), without shifting the transition barycenter. 
The total accuracy of our measurements, summarized in 
TableQJ ranges from lxlO~ n to 5xl0" 12 , which is cur- 
rently the best reported result for any optical transition 
in 3 He. 

An independent check for the accuracy of our mea- 
surements is made by extracting the known value of 
the hyperfine splitting (HFS) of the metastable 2 3 5 
state from the transitions in Table |TJ The two dif- 



ferences P, 



,1/2 



and Py^-Pf 7 , 2 ,, , yield the val- 



p3/2 , „l/2 _p3/2 

0.1/2 ^0.1/2 ailCL 1,1/2 "1,1/2 

ues of 6 739 701.5 (3.0) and 6 739 701.3(3.1) kHz, re- 
spectively, which are consistent with each other and 
in perfect agreement with the more accurate result of 
6739701.177(16) kHz [l(J. 

Comparison of our measurements with the previous ex- 
perimental results is given in Table HI The centroid val- 
ues of the 2 3 P and 2 3 S' energies are defined as an average 
over all fine and hyperfine sublevels, 



E(2 3 L) = 



J: jf (2F+1)E(2 3 Lj, f ) 



(1) 



(27+l)(25+l)(2L + l) ' 
where 2 3 P = 2 3 6* and 2 3 P for L= and 1, respectively. 
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FIG. 2. (color online) Day-by-day correction of the 

1/2 

RS+LS+2DS shift for the measurements of the f 2 3/2 
perfine transition. Left graph: the squares (black) are the 
measured data, and the circles (red) are the corrected data. A 
different mean frequency was subtracted to bring data in the 
same frequency vertical scale for a clear comparison. Right 
graph: bar graph distribution of the measured (dashed black 
bars) and corrected (fill red bars) frequencies. 



TABLE II. Comparison with prior measurements, in kHz. 
3 He 2 3 P-2 3 S' centroid 

4 He 2 3 P-2 3 S' centroid 



4 He P 2 - 3 He F^J? 



0,1/2 



3He <i 2 /2- 4He ^1 



276 702 827 204.8 (2.4) 
276 702 827 145 (77)^ 

276 736 495 649.5 (2.1) [7] 
276 736 495 580(70) [11] 

810 594.3 (3.3^ 
810 599(3) [3J 
810 608(30) [12] 

1480 582.1 (3.2) 

1 480 573 (30) [12] 



Evaluated by combining the 4 He P 2 frequency [Tl| . the 
4 He- 3 He Pl~P^l j2 interval 0, the 3 He 2 3 S HFS 0, and the 
2 3 PHFS HI. 

Evaluated by combining the fj^jyj i n t erv al from Table[T] the 
P Q frequency 0, and the 2 3 P -2 3 P 2 interval 0. 



To the best of our knowledge, there are no published 
measurements of 2 3 P-2 3 5 HFS frequencies. Therefore, 
we obtain the "previous experimental value" of the 3 He 
2 3 P-2 3 5' centroid frequency as a combination of several 
experiments and the calculated HFS intervals (see foot- 
note inTableHH. The previous result is in agreement with 
our measurement but 33 times less accurate. The 4 He 
2 3 P-2 3 5' centroid was measured by us previously 0], in 
agreement with the independent determination by Shiner 
et at 11]. In order to check the consistency of our 
present measurements on 3 He with our previous mea- 



surements for 4 Hc 0], in Table HI a comparison of the 
frequency differences of the 3 He P^^m an d the 4 Hc Pi 



TABLE III. 4 He- 3 He isotope shift of the 
for the pointlike nucleus, in kHz. m r is the 
M is the nuclear mass. 



Contribution 



2 3 P-2 3 S 



centroid energies, 
reduced mass and 

2 1 S-2 3 S 



m r a 2 (m r /M) 
m r a 2 (m r /M) 2 
m r a 2 (m r /M) 3 
m r a 4 

m r a 4 (m r /M) 

m r a 4 (m r /M) 2 

ma 5 (m/M) 

ma 6 (m/M) 

Nuclear polarizability 

HFS mixing 

Total 

Other theory [15, l&£ 



12412 458.1 
21243 041.3 
13 874.6 
4.6 
17 872.8 
-20 082.4 
-3.0 
-60.7 
-15.5(3.9) 
-1.1(1) 
54.6 

33 667143.2 (3.9) 
33 667146.2 (7) 



8 632 567.86 
-608 175.58 

7 319.80 
-0.30 

8 954.22 
-6 458.23 

-1.84 
-56.61 
-2.75 (69) 
-0.20 (2) 
-80.69 
8 034 065.69 (69) 
8 034 067.8(1.1) 



a Corrected by adding the triplet-singlet HFS mixing. 

and Pi intervals with independent measurements 0, [l2| 
is reported. 

The difference of our results for the 2 3 P-2 3 S I centroid 
energy in 3 He and 4 He yields the experimental value of 
the isotope shift (IS). Combined with theoretical calcu- 
lations, the experimental IS can be used 14] to deter- 
mine the difference of the squared nuclear charge radii, 
6r 2 



„2{3 



He) -r 2 ( 4 He). 
Numerical results of our calculation of the IS of the 
2 3 P-2 3 5'and 2 1 S- 2 3 S transitions for the point nucleus are 
presented in Table HTT1 As compared to the previous eval- 
uations 15, HH, the higher-order recoil [m r a 2 (m r /M) 3 
and m r a 4 (m r /M) 2 ] and the nuclear polarizability cor- 
rections were accounted for and the higher-order QED 
effects [ma 6 (m/M)] were estimated more carefully. The 
calculation extends our previous works 17 , 18{ ; its details 
will be reported elsewhere. The total uncertainty comes 
from the uncalculated higher-order QED effects. Note 
that it is much larger than the one reported previously 
in Ref. 0. 

Our definition of the IS differs from that used previ- 
ously 0, [If by the fact that we average out not only 
the hyperfme but also the fine-structure splitting. The 
advantage of using the centroid energy is that theory be- 
comes much more transparent and can be directly com- 
pared to the experiment. 

The difference of the theoretical point-nucleus results 
in Table ITTll and the experimental IS comes from the finite 
nuclear size (FNS) effect, which can be parametrized as 



1 - (Zaf \n(Zar) + (Za) 2 f le i 



(2) 



where / re i is the relativistic correction beyond the lead- 
ing logarithm. The FNS contribution can be represented 



4 



as J-Efns = C r 2 , where the coefficient C, according to 
the above equation, depends very weakly on r. Our cal- 
culated values for the coefficient C are 



particular importance, as the next goal of the muonic 



C(2 3 P - 2 3 S) = -1212.2 (1) kHz/fm 2 , 
C(2 1 S - 2 3 S) = -214.69 (1) kHz/fm 2 , 



(3) 
(4) 



which can be compared with the previous results C(2 3 P— 
2 3 S) = -1209.8 pl] and C(2 1 S - 2 3 S) = -214.40 [U- 

At present, there are three independent measurements 
of the 4 He- 3 He IS that can be used to infer 5r 2 with 
a comparable accuracy, our experiment and those of 
Refs. [2, 3]. Our theory, summarized in Table Hill allows 
us to extract the charge radius difference 5r 2 consistently 
from the present experiment and that of Ref . |2J . The re- 
sults are 

5r 2 (this work) = 1.074 (3) fm 2 , (5) 
<5r 2 ([2], reevaluated) = 1.028 (11) fm 2 . (6) 

The 5r 2 value of Eq. (j6|) is by about la larger than that 
given in Ref. 0], 1.019 (11) fm 2 , because of the change in 
the theoretical IS value. Using Eq. §5§ and the nuclear 



charge radius of 4 He 2l| , we obtain the root-mean-square 
radius of the 3 He isotope to be 1.975 (4) fm. 

The results of Eqs. (J5]) and dU can be also compared 
with the determination by Shiner et al., 



fr 2 (0) = 1-059 (3)fm 2 



(7) 



which used the older isotope shift theory. We do not 
reevaluate this result here, since it would require im- 
provement in theoretical predictions for HFS intervals. 

We observe, as shown in Fig|3l that the above results 
for the radius difference Sr 2 are inconsistent with each 
other. In particular, a 4 a discrepancy is present between 
our value and that of Ref. 0] , for which we do not have 
a satisfactory explanation at present. Note that both ex- 
periments use the same OFCS assisted laser technology, 
that most of the theoretical contributions are checked by 
independent calculations, and that the determination of 
the charge radius difference is now performed consistently 
within the same theory. The observed discrepancy may 
be in principle explained by some hidden systematics in 
experiments or by yet unknown effects in the electron- 
nucleus interaction. 

The possibility that some additional effects beyond the 
standard QED exist has been discussed in the literature 
ever since the muonic hydrogen experiment [H raised 
what is now known as the proton charge radius puzzle. 
One of the ways for solving this puzzle is to investigate 
similar systems, aiming to confirm or to disprove the dis- 
agreement observed for hydrogen. In the present Letter, 
we report a 4cr discrepancy for the nuclear charge radius 
difference of 3 He and 4 He. 

Precision spectroscopic determination of the nuclear 
charge radii of the helium isotopes becomes today of 



■» From 2'P-2'S IS [this work] 
-From 2'S-2 3 S IS [2], reevaluated as Eq.(6) 



From electron-scattering [22,23] 



From nuclear theory [20,21]'- 



Sr 2 [ 3 He- 4 He] (fm 2 ) 

FIG. 3. (color online) Different determinations of the differ- 
ence of the squared nuclear charge radii for 3 He and 4 He. 



hydrogen group from the Paul Scherrer Institute is the 
measurement of the Lamb shift in muonic helium 



24] 



This experiment will provide an independent determina- 
tion of the charge radii of helium isotopes and will allow 
one to compare the results obtained by the spectroscopy 
of the electronic and muonic atoms, thus hopefully shed- 
ding light on the proton charge radius problem and on 
the discrepancy for the helium charge radius difference. 
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